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We consider axially symmetric, rotating boson stars. Tlieir fiat space limits represent spinning 
Q-balls. We discuss their properties and determine their domain of existence. Q-balls and boson 
stars are stationary solutions and exist only in a limited frequency range. The coupling to gravity 
gives rise to a spiral-like frequency dependence of the boson stars. We address the flat space limit 
and the limit of strong gravitational coupling. For comparison we also determine the properties of 
spherically symmetric Q-balls and boson stars. 
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^ ; I. INTRODUCTION 

Non-topological solitons [1] or Q-balls [2] represent stationary localized solutions in flat space possessing a finite 
■ mass. In the simplest case they arise, when a complex scalar field has a suitable self-interaction, mimicking the 
\ interaction with other fields [1] . The global phase invariance of the scalar field theory is associated with a conserved 

^\ • charge Q, corresponding for instance to particle number [1]. Q-balls thus represent non-topological solitons with 

PsJ ' charge Q. 

Spherically symmetric Q-balls exist only in a certain frequency range, Wj^in < i-^s < ^^maxi determined by the 
properties of the potential [1-3] . At a critical value of the frequency, both mass and charge of the Q-balls assume 
their minimal value, from where they monotonically rise towards both limiting values of the frequency. Considering 
the mass of the Q-balls as a function of the charge, there are thus two branches of Q-balls, merging and ending at 
the minimal charge and mass. Q-balls are stable along the lower branch, when their mass is smaller than the mass of 
• Q free bosons [1]. 

' When gravity is coupled to Q-balls, boson stars arise [3, 6-8]. The presence of gravity has crucial influence on 
lO . the domain of existence of the classical solutions. Stationary spherically symmetric boson stars also exist only in 
' a limited frequency range, Wo(k) < uJs < Wmaxj where n denotes the strength of the gravitational coupling. They 
, show, however, a different type of critical behaviour. For the smaller frequencies the boson stars exhibit a spiral-like 
O^' frequency dependence of the charge and the mass, approaching finite limiting values at the centers of the spirals. 
^ When the maximal value of the frequency is approached, the charge and the mass of the boson stars tend to zero 
[3, 6]. 

^ I We here determine the dependence of the stationary spherically symmetric boson stars on the strength of the 
gravitational coupling k and, in particular, determine the domain of existence of the boson stars as a function of k, 
including the limit k oo and the flat space limit. 

Recently, the existence of rotating Q-balls was demonstrated [4] . These correspond to stationary localized solutions 
in flat space possessing a finite mass and a finite angular momentum. Interestingly, their angular momentum is 
quantized, J = nQ [4, 5]. Possessing even or odd parity, their energy density forms one or more tori. 

Rotating boson stars exhibit angular momentum quantization as well, J = nQ [5]. Previously they have been 
obtained only for potentials without solitonic flat space solutions, and only for a restricted frequency range [5, 9- 
11]. Here we consider rotating boson stars which possess flat space counterparts. Our main objectives are to clarify, 
whether rotating boson stars exhibit an analogous frequency dependence as non-rotating boson stars, and to determine 
their domain of existence. For that purpose, we compute numerically sequences of rotating boson stars for constant 
values of the gravitational coupling strength, focussing on fundamental solutions with rotational quantum number 
n = 1 and even parity. 

In section II we recall the action, the general equations of motion and the global charges. In section III we present 
the stationary axially symmetric ansatz for the metric and the scalar field, we evaluate the global charges within 
this Ansatz, and present the boundary conditions for the metric and scalar field functions. We discuss stationary 
spherically symmetric Q-ball and boson star solutions in section IV, and present our results for rotating Q-ball and 
boson star solutions in section V. Section VI gives our conclusions. In the Appendices we address two auxiliary 
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functions employed in the numerical integration (A), we present the systems of differential equations (B) and the 
components of the stress-energy tensor (C). 



II. ACTION, EQUATIONS AND GLOBAL CHARGES 
A. Action 



We consider the action of a self-interacting complex scalar field $ coupled to Einstein gravity 



S = 



I 



gd X , 



ISttG 2' 

where R is the curvature scalar, G is Newton's constant, the asterisk denotes complex conjugation, 



and U denotes the potential 



U{\^\) = A|$|2 - a|$|2 + {,) = A - a<f>^ + bc/)'^) 



(1) 



(2) 



(3) 



with |4>| = (j). The potential is chosen such that nontopological soliton solutions [1], also referred to as Q-balls [2], 
exist in the absence of gravity. As seen in Fig. 1, the self- interaction of the scalar field has an attractive component, 
and the potential has a minimum, C/(0) = 0, at i> = and a second minimum at some finite value of |$|. The boson 
mass is given by me — -\/A6. 
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FIG. 1: The potential U{4>) is shown for A = 1, a = 2 and b = 1.1 resp. 6 = 1. 



B. Equations 



Variation of the action with respect to the metric leads to the Einstein equations 



with K = SttG and stress-energy tensor T^v 



Tfxv — Qui/Lm — 2 



dL 



(4) 

(5) 
(6) 
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Variation with respect to the scalar field leads to the matter equation, 

where □ represents the covariant d'Alembert operator. Equations (4) and (7) represent the general set of Einstein- 
Klein-Gordon equations. 



C. Global Charges 

The mass M and the angular momentum J of stationary asymptotically flat space-times can be obtained from their 
respective Komar expressions [12], 



M 

and 



J=-Z^ I R^..n''r]''dV . (9) 



SttG 



E 



Here S denotes an asymptotically flat spacelike hypersurface, is normal to S with n^n^ = — 1, dV is the natural 
volume element on E, ^ denotes an asymptotically timelike Killing vector field and r] an asymptotically spacelike 
Killing vector field [12]. Replacing the Ricci tensor via the Einstein equations by the stress-energy tensor yields 

M = 2 ^ (^T^, - i g^, T^-y^ n^'CdV , (10) 

and 

J = -J^(t^.-1 ff^. v) ^^^'^'^v ■ (11) 

A conserved charge Q is associated with the complex scalar field since the Lagrange density is invariant under 
the global phase transformation 

$ ^ $e*" (12) 

leading to the conserved current 

j** = _j ($*a''$ - ^a****) , i% = 0. (13) 

III. ANSATZ AND BOUNDARY CONDITIONS 
A. Ansatz 

To obtain stationary axially symmetric solutions, we impose on the space-time the presence of two commuting 
Killing vector fields, ^ and ij, where 

^ = dt, r] = d^ (14) 

in a system of adapted coordinates {t, r, 9, ip}. In these coordinates the metric is independent of t and (p, and can be 
expressed in isotropic coordinates in the Lewis-Papapetrou form [13] 



g {dr'^ + dG'^) + sin^ 6 (dtp-'^ dt) ' 
The four metric functions /, I, g and w are functions of the variables r and 9 only. 



(15) 
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The symmetry axis of the spacetime, where r] = 0, corresponds to the 2;-axis. The elementary flatness condition 
[14] 

^^ = 1, X = rj^rj, (16) 
then imposes on the symmetry axis the condition [13] 

g\0=o = g\0=-K = 1 • (17) 

For the scalar field $ we adopt the stationary ansatz [5] 

^{t,r,e,ip) = (t>{r,e)e^'*+''"^ , (18) 
where (j){r, 6) is a real function, and and n are real constants. Single-valuedness of the scalar field requires 

$(99) = $(27r + ip) , (19) 

thus the constant n must be an integer, i.e., n = 0, ±1, ±2, .... We refer to n as rotational quantum number. When 
n ^ 0, the phase factor exp (imp) prevents spherical symmetry of the scalar field 

Thus to obtain stationary axially symmetric boson stars a system of five coupled partial differential equations needs 
to be solved. This set of equations is presented in Appendix B 2. In contrast, for the Q-balls of flat space the metric is 
the Minkowski metric, i.e.,/ = I = g = 1, uj = 0. Here, at least in principle, only a single partial differential equation 
for the scalar field function needs to be solved. 

For stationary spherically symmetric boson stars n = 0, and the scalar field function depends only on the radial 
coordinate, ([> = 4>{r). The metric then simplifies as well, since g = 1 and lu = 0, and the non-trivial metric functions 
/ and I depend only on the radial coordinate, f = f{r), I = l{r). Thus for stationary spherically symmetric solutions 
one obtains a much simpler system of three coupled ordinary differential equations. This set of equations is presented 
in Appendix B 1. 



B. Mass, angular momentum and charge 

The mass M and the angular momentum J can be read off the asymptotic expansion of the metric functions / and 
CO, respectively, [15] 



/ = 1 



2MG 



2JG 



(20) 



I.e., 



M = ^ lira r^dr f , J = lim r^u . 



(21) 



This is seen by considering the Komar expressions Eqs. (8) and (9), with unit vector n'* = (1, 0, 0, a;/r)/v7) and 
volume element dV = Ij y/] \g\^l'^ dr dB d(p, leading to [16] 



M 



Rly/^drdOdip 



lim — 

r->oo 8nG 



f 

Jo 



f 



df I 



( duj 



—r sin0 I — -sin Ouj ( — 



dr f 



\dr r 



de 



(22) 



and similarly 



J = lim 



27r 



00 167rG 



Jo 



-j^r sm 0\^^-r— 



de 



(23) 



Insertion of the asymptotic expansions of the metric functions then yields expressions (21). 
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Alternatively, the mass M and the angular momentum J can be obtained by direct integration of the expressions 
(10) and (11), 

M= Jj2T;-S;T:j)n.^'^dV , 

= J (2T/-T;) \g\'/^drd0d^ , (24) 

corresponding to the Tolman mass, and 

J = - j Tl\g\^/^drded<fi . (25) 
The conserved scalar charge Q is obtained from the time-component of the current, 

Q = - j f \g\^/^ drdddifi 

Prom Eq. (25) for the angular momentum J and Eq. (26) for the scalar charge Q, one obtains the important 
quantization relation for the angular momentum, 

J = nQ, (27) 

(Icrivcd first by Schunck and Miclkc [5], by taking into accomit that T^^ = nj*, since = m$. Thus a spherically 
symmetric boson star has angular momentum J = 0, because n = 0. 



C. Boundary conditions 

The choice of appropriate boundary conditions must guarantee that the boson star solutions are globally regular 
and asymptotically flat, and that they possesses finite energy and finite energy density. 

For spherically symmetric boson stars boundary conditions must be specified for the metric functions /(r) and l{r) 
and the scalar field function (f){r) at the origin and at infinity. At the origin one finds the boundary conditions 

drf\r=0 = , drl\r=0 = , 9r0|r=O = . (28) 

Note, that for spherically symmetric boson stars the scalar field has a finite value (po at the origin, 

<P{r) =<Po+ O(r') . (29) 

For r — > oo the metric approaches the Minkowski metric rjap and the scalar field assumes its vacuum value $ = 0. 
Accordingly, we impose at infinity the boundary conditions 

/Ir^oo = 1 , ^Ir^oc. = 1 , (p\r^oo=0. (30) 

For rotating axially symmetric boson stars appropriate boundary conditions must be specified for the metric func- 
tions f{r,9), l{r,9), g{r,6) and uj{r,9), and the scalar field function (l>{r,9) at the origin, at infinity, on the positive 
2;-axis (^ = 0), and, exploiting the reflection symmetry w.r.t. 9 ^ w — 6, in the xy-plane {9 = 7r/2). At the origin we 
require 

drf\r=0 = , drl\r=0 = , g\r=0 = 1 , Oj\r=0 = , 0|r=O = . (31) 

At infinity the boundary conditions are 

/Ir^oo — 1 , ^Ir^oo — 1 , g\r^oo — 1 , ^|r^oo — , (j)\r^oo = , (^^) 

and for 9 = and 9 = 7r/2, respectively, we require the boundary conditions 

90/|e=o = , del\e=o = 0, g\0=o = 1 , i9ew|e=o = , (l)\e=o = 0, (33) 
and for even parity solutions 

9e/|e=77/2 = , d0l\g^^/2 = O, dgg\0=T,/2 = , 9ew|e=7r/2 = , de(p\0=T,/2 = , (34) 
while for odd parity solutions <P\0=t!-/2 = 0. 
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IV. STATIONARY SPHERICALLY SYMMETRIC SOLUTIONS 

Stationary spherically symmetric solutions are obtained, when n = 0. The set of coupled non-linear ordinary 

differential equations, given in Appendix Bl, is solved numerically [17], subject to the above boimdary conditions, 
Eqs. (28)-(30). Because of the power law fall-off of the metric functions, we compactify space by introducing the 
compactified radial coordinate 

f = . (35) 

1 + r ^ ^ 

The numerical calculations employ a collocation method for boundary-value ordinary differential equations, using the 
damped Newton scheme for a sequence of meshes, until the required accuracy is reached [17]. 

A. Solutions in flat space: Q-balls 

Spherically symmetric Q-balls have been studied before (see e.g. [1, 2, 4], or [3] for a review). We here review their 
main features, to be able to better demonstrate the effects of gravity and rotation. 
Recalling first the mass M of the Q-balls, 

/.oo poo 

M = 4n TttT^dr = An [u"^^ cp"^ + (j)'^ + U{(p)] dr , (36) 
Jo Jo 

where the prime denotes differentiation with respect to r, and and their charge Q, 

Q{uJs) = Sttojs (fi^r'^dr , (37) 
Jo 

we follow the discussion in [4] to obtain the limits for the frequency 

while correcting for the missing factor of two in the field equation [18, 19]. The equation of motion for the scalar field 
[18], 

= r + -<^'-'/-^+-^.^, (39) 
r 2 d(p 



is equivalent to [18] 

! / —dr, (40) 
Jo r 



111 rr 



where £ is an integration constant, and effectively describes a particle moving with friction in the potential [18] 

V{<l>) = \^lcf'-\u{cj^) . (41) 
The first necessary condition for Q-balls to exist is then V" (0) < [4] and yields the maximal frequency Wmax [18] 

< = \u"{0) = \h, (42) 

while the second condition is that V{(1)) should become positive for some non-zero value of (j) [4], i.e., [18] 

ujI > Lol^ = mm [i7(<^)/</.2] = A (^6 - ^) . (43) 

Turning now to the Q-ball solutions, we specify the potential parameters as [4] 

A = 1 , a = 2 , 6 = 1.1 . (44) 
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FIG. 2: Left: The charge Q is shown as a function of the frequency lVs for fundamental Q-balls (fc = 0) and their first radial 
excitations {k = 1). Also shown are the limiting values of the frequency, a;niin and Wmax- Right: The mass M is shown as a 

function of the charge Q for fundamental Q-balls (fc — 0) and their first radial excitations (fc = 1). Also shown is the mass 
for Q free bosons, M = tubQ- The upper branches of the mass M are not discernible (on this scale) from the mass of Q free 
bosons. 



Fixing the value of uJs in the allowed range, one obtains a sequence of Q-ball solutions, consisting of the fundamental 
Q-ball and its radial excitations [4]. The boson function 4> of the fundamental Q-ball has no nodes, while it has k 
nodes for the fc-th radial excitation. 

Focussing on the fundamental Q-ball solutions and their first radial excitations, we exhibit in Fig. 2 the dependence 
of the charge Q on the frequency tOg- As seen in the figure, at a critical value Wcr the charge assumes its minimal 
value Qcr- The charge diverges both when lOs i^'min find when lOs — > Wmax [!]■ 

Fig. 2 also exhibits the mass M as a function of the charge Q. The region close to the critical value Qcr is exhibited 
in Fig. 3 for the fundamental Q-balls and their first radial excitations. The lower branches correspond to values of 
the frequency LOg < Ucr, while the upper branches represent the values lUs > (^cr- When the mass is smaller than the 
mass of Q free bosons, M < ttibQ, the solutions are stable [1]. When iVg Wmax, the upper branches approach the 
mass M = itibQ of Q free bosons from above [1]. 
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FIG. 3: The mass M is shown as a function of the charge Q for fundamental Q-balls (fc = 0, left) and their first radial 
excitations (fe = 1, right) in the region close to their corresponding critical values Qcr. Also shown is the mass for Q free 
bosons, M = tubQ- 



The scalar field function (p and the energy density of the fundamental Q-balls are shown in Fig. 4 for several 
values of the frequency iVg including Wcr- The energy density of the Q-balls is shell-like. Along the upper branch, with 
increasing charge and mass the maximum of the energy density decreases while moving outwards, whereas along the 
lower branch, with increasing charge and mass the maximum of the energy density increases while moving outwards, 
being strongly correlated with the steep fall-off of the scalar field function. Radially excited Q-balls with k nodes 
possess k + 1 energy shells [4] . 



8 




FIG. 4: The scalar field </)(r) and the energy density Tu(r) are shown for fundamental Q-balls (fc = 0) for several values of the 
frequency ojs including Wcr- For uig = 0.6 and u)s = 0.4152 the scaled energy density Ttt{r)/4 resp. Ttt{r) /25 is shown. 

B. Solutions in curved space: boson stars 

When the scalar field is coupled to gravity, boson stars arise (see e.g. [7, 9] for reviews). Spherically symmetric 
boson stars, based on a self- interacting boson field with a potential U, Eq. (3), have been considered by Friedberg, 
Lee, and Pang [6]. For their choice of the potential parameters the minima of the potential are degenerate, yielding 
for the minimum value of the frequency uj^in = 0. In contrast, for our choice of parameters, Eq. (44), the potential 
has a global and a local minimum, and Wmin > 0. The different form of the potential has consequences for some 
features of the boson stars, as discussed below. 

To demonstrate the effects of gravity on the spherically symmetric solutions, we exhibit in Fig. 5 the charge Q as 
a function of the frequency for the fundamental boson star solutions at a given value of the gravitational coupling 
constant k. 




FIG. 5: The charge Q is shown as a function of the frequency uJs for fundamental boson stars (fe = 0) in the full range of 
existence (left), and in the frequency range of the spiral (right), for the gravitational coupling constant k = 0.2. 

For solutions in curved space the frequency is also bounded from above by Wmax, Eq. (42), since the scalar 
field exhibits asymptotically an exponential fall-off only for oJs < Wmax- However, for the smaller values of ojg a new 
phenomenon occurs, as compared to flat space [6], where the solutions approach monotonically the limiting lower 
value Wmin, Eq. (43). Denoting the minimal value of the frequency for boson stars by wo('«), we observe that it differs 
from Wmin, Bxcept for a single value of the gravitational coupling k. For larger values of k, u}o{i^) > ^^min, while for 
smaller values of k, wo(k) < ujmin- (Note, that the latter case docs not occur for a potential with degenerate minima, 
since there Wmin = 0.) Moreover, in the presence of gravity the solutions do not approach monotonically the minimal 
value wo{k,). Instead one observes an inspiralling of the boson star solutions towards a limiting solution at the center 
of the spiral at a frequency Lj\i^{K) > ujo{k) [6]. 

In Fig. 6 we exhibit the charge Q as a function of the frequency oJs > for the fundamental boson star solution (fc = 0) 
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for several values of the gravitational coupling constant k, and also for the first radial excitation (fc = 1) for k = 0.2. 
While the presence of a spiral is a genuine property of boson stars, the location and the size of the spiral depend on 
the gravitational coupling strength k and on the node number k. 
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FIG. 6: The charge Q is shown as a function of the frequency lOs for fundamental boson stars {k = 0) for the values of the 
gravitational coupling constant k = 0.1, 0.2, 0.4, 1.0 (left), and for their first radial excitations (fc = 1) for k = 0.2 (right). 



The mass M has an analogous dependence on the frequency ojs as the charge Q, as seen in Fig. 7, where we also 
exhibit the frequency dependence of the value of the; scalar field at the origin 4>{0). The endpoints of 4>{0) in the 
figure correspond to the numerical values obtained closest to the centers of the spirals. We cannot decide numerically, 
however, whether the limiting solutions have indeed finite values ^iim(O) [20]. 
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FIG. 7: The mass M (left) and the value of the scalar field at the origin ^(O) (right) are shown as functions of the frequency LOa 
for fundamental boson stars {k = 0) for the values of the gravitational coupling constant k = 0.1, 0.2, 0.4, 1 (left) resp. k = 0.02, 
0.1, 0.2, 0.4, 1 (right). 



When the mass M is considered as a function of the charge Q, we observe a cusp structure [6], as illustrated in 
Fig. 8. For smaller values of k (e.g., k = 0.2) we observe two sets of cusps. The first set of cusps is related to the 
single cusp present in flat space, which occurs at the minimal value of the charge Qmin, where the upper and lower 
branch of the Q-ball merge. But whereas the upper branch extends infinitely far in flat space, it extends in curved 
space only up to a second cusp, reached at a second critical value of the charge, where this upper branch merges with 
a third branch extending down to zero. For larger values of k (e.g., k = 1) this set of two cusps, being a remainder 
of the flat space solutions, is no longer present. 

The second set of cusps is related to the spirals, present only in curved space. Labelling these cusps of the spirals 
consecutively A'' = 1, 2, we exhibit the charge Q, the mass M, and the frequency Us of the boson star solutions at 
these cusps in Table I, for several values of the gravitational coupling constant k. 

Focussing on the limiting solutions at the centers of the spirals, we exhibit in Table II the limiting values of the 
frequency wum, the charge Qnm and the mass Mum for several values of the gravitational coupling constant k. For 
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FIG. 8: The mass M is shown as a function of the charge Q for fundamental boson stars (k = 0) for the values of the 
gravitational coupling constant k = 0.2 (left) and k = 1 (right). 



N 


Q-v 


M 






N 


Qn 


M 




1 


685.5000 


351.6320 


0.34338 




1 


119.2400 


122.9620 


0.89402 


2 


331.3000 


221.3990 


0.42424 




2 


78.1580 


88.3579 


0.89343 


3 


391.9200 


246.5870 


0.40352 




3 


84.8980 


94.4262 


0.89000 


4 


378.2500 


241.0480 


0.40776 




4 


83.4960 


93.1796 


0.89067 


5 


380.8700 


242.1160 


0.40693 




5 


83.7470 


93.4033 


0.89047 


6 


380.4000 


241.9250 


0.40708 




6 


83.7090 


93.3694 


0.89048 


N 


Qn 


M 


UJs 




N 


Qn 


M 


UJs 


1 


29.3040 


28.1273 


0.63000 




1 


12.2180 


12.5360 


0.93000 


2 


15.7880 


19.1084 


0.76800 




2 


5.6059 


7.24552 


0.86000 


3 


17.8630 


20.6836 


0.74018 




3 


6.3314 


7.87359 


0.85305 


4 


17.4580 


20.3833 


0.74497 




4 


6.1911 


7.75446 


0.85290 


5 


17.5280 


20.4361 


0.74416 




5 


6.2148 


7.77456 


0.85304 


6 


17.5170 


20.4273 


0.74429 




6 


6.2108 


7.77104 


0.85302 



TABLE I: Physical characteristics of the first six cusps of the fundamental boson star solutions {k = 0) for gravitational coupling 
constants k = 0.1 (upper left), k = 0.2 (upper right), k = 0.4 (lower left), k = 1 (lower right). 



small gravitational coupling we extract the following K-dependence, 

Qiim«'/' = c;? + c?«:i/2+c?^ + 0(^3/2) ^ 

Mh^,«3/2 ^ ^ ^M^l/2 ^ ^M^ ^ o(^3/2) ^ (45) 

illustrated in Fig. 9. For small gravitational coupling wc thus obtain a different K-dcpendcnce, as compared to [6]. In 
particular, in the limit k — > the limiting frequency assumes a finite value, < wiini(O) < Wmin- We attribute 
this difference to the fact, that for a potential with degenerate minima, as employed in [6], Wmin = 0. 

To address the K-dcpendcncc of the domain of existence of the (fundamental) boson star solutions, corresponding 
to the interval [[jJo(k), tiJniax], we exhibit in Fig. 10 the charge Q as a function of the frequency lUs for a large set of 
values of the gravitational coupling constant k, including the flat space limit. Interestingly, for values of the frequency 
UJs close to Wmax, the flat space values are not approached monotonically from below with decreasing k. Here a weak 
coupling to gravity can lead to an increase of the mass of the boson star solutions. 

As K increases, cjo(k) tends to a finite value ujq^oo), smaller than Wmax- Fixing ujg G [liJo(oo), o^max] wc observe that 
the scalar field scales like 1/-\/k for large values of k. In order to obtain the solutions in the limit k ^ 00, we therefore 
introduce the scaled scalar field (j){r) = s/K(j){r). Substituting (j){r) in the field equation and taking the limit k ^ 00, 
we find that all terms non-linear in (f){r) vanish. Similarly, in the Einstein equations all terms of higher than second 
order in (j}{r) vanish and the dependence on k cancels. Thus we end up with a set of differential equations identical 
to the original one, except that k = 1 and U{(j)) = Xb(j)^ . Solving this set of differential equations for fixed Xb and 
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Qlim 




1.0000 


0.85302 


6.2114 


7.7718 


0.4000 


0.74427 


17.5180 


20.4285 


0.2000 


0.56135 


70.8330 


63.2038 


0.1000 


0.40707 


380.4600 


241.9480 


0.0200 


0.25981 


1.136 • 10^ 


4433.790 


0.0100 


0.23529 


4.011 • 10* 


1.401 • 10* 


0.0020 


0.20754 


5.954 ■ lO'-' 


1.802 ■ lO"^ 


0.0010 


0.20172 


1.797 ■ lO'' 


5.262 ■ 10® 


0.0002 


0.19451 


2.181 • lO'' 


6.119 • 10® 


0.2 • 10"® 


0.18924 


2.326 • 10^° 


6.316 • 10^ 



TABLE II: Physical characteristics of the hmiting solutions at the centers of the spirals for the fundamental boson stars (k = 0) 
for a set of decreasing values of the gravitational coupling constant k. 
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FIG. 9: The dependence of the frequency wum, the mass Mum and the charge Qum of the limiting solutions at the centers of 
the spirals on the gravitational coupling constant k for fundamental boson stars (fc = 0) for small values of k. 
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FIG. 10: The charge Q (left) and the mass M (right) are shown as functions of the frequency uig for fundamental boson stars 
(fc = 0) for a set of decreasing values of the gravitational coupling constant k. Also shown are the limiting flat space values. 



varying Ug, we again find a spiral pattern for the (scaled) charge Q = kQ, shown in Fig. 11. This analysis yields 

a;o(oo) = 0.805 , a;iim(oo) = 0.883 . 



Considering now the domain of existence in the limit k ^ 0, we note that we must distinguish two intervals here, 
[u)o{k) , u)^i-a\ and [wminj'^inax]- Wc observe from Fig. 10, that the limiting flat space values of the charge Q are 

approached in a continuously increasing interval, extending to the full interval [t^'min, ^max] in the limit n ^ 0. In the 
interval [a;o(0), Wmin] there are no flat space solutions, possessing flnite charge or mass. As noted above, the limiting 
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FIG. 11: Left: The scaled charge kQ is shown as a function of the frequency uJs for fundamental boson stars (fe = 0) in the limit 
K ^ oo. For comparison the scaled charge kQ is also shown for k = 1. resp. k —* (right). Right: The scaled charge k^^^Q 
is shown as a function of the frequency oj, for fundamental boson stars (A; = 0) for k = 0.002, 0.0002, 0.00002, and 0.000002. 
The asterisk marks Wmin- 



values of the charge Qnm and the mass Miim at the centers of the spirals diverge with when k ^ 0. We therefore 

determine the properties of the sohitions in the limit k ^ in the interval [cl)o(O), tiJmin] by scaling the charge and the 
mass with k^/^ for a sequence of solutions, corresponding to decreasing values of k. Convergence towards the limiting 
values of the scaled charge is demonstrated in Fig. 11 [21]. Thus also in the limit k; — > we find a spiral pattern for 
the (scaled) charge Q = k^^^Q. This analysis yields 

a;o(0) = 0.172 , a;iim(0) = 0.189 . 

We remark that for a potential with degenerate minima, as employed in [6], only the single interval [wmin, Wmax] needs 
to be considered, since Wmin = 0, implying apparently a;o(0) = a;iiin(0) = Wmin = 0. 



V. ROTATING AXIALLY SYMMETRIC SOLUTIONS 



Rotating axially symmetric solutions are obtained, when n ^ 0. We solve the set of coupled non-linear elliptic 

partial differential equations, given in Appendix B2, numerically [22], subject to the above boimdary conditions, 
Eqs. (31)-(34), employing the compactificd radial coordinate, Eq. (35). The numerical calculations are based on the 
Newton-Raphson method. The equations are discrctizcd on a non-equidistant grid in f and 6. Typical grids used 
have sizes 100 x 20, covering the integration region < f < 1 and < ^ < 7r/2. 



A. Rotating Q-balls 

The existence of rotating Q-balls has been shown by Volkov and Wohnert [4]. Based on the ansatz (18) for the 
scalar field $ [5], one obtains for rotating Q-ball solutions the field equation [18] 

2d_ J_d^ COS0 d A rh - l^EMl 

the mass [4] 

M = 27r r r (lo^,4>^ + {dr<pf + \{de<Pf + , + C/(0)) dr smOdO , (47) 

Jo Jo \ ^ Hsm^6' J 

and the charge 

Q = 4nuJs / (j?r^dr sme dO . (48) 
Jo Jo 

Their angular momentum satisfies the quantization relation J = nQ. 
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In their pioneering study [4] Volkov and Wohncrt show, that for a given value of n there arc two types of solutions, 
possessing different parity. They exhibit examples of fundamental rotating Q-balls with quantum numbers n = 1 — 3, 
and both even and odd parity. They do not study the frequency dependence of rotating Q-ball solutions, however. 

We illustrate the scalar field 0(r, ^) and the energy density Ttt{r,0) of fundamental rotating Q-balls with charge 
Q = 410, quantum numbers n = 1 and 2, and even parity in Figs. 12 and 13. For rotating Q-balls the scalar field (f) 
must vanish at the origin, as seen in Eq. (47). The energy density Ttt of even parity rotating Q-balls is torus-like. 
(For odd parity a double torus arises [4].) The mass M, the angular momentum J and the frequency Wg of these 




FIG. 12: The scalar field 4>{r, 9) for fundamental rotating Q-balls with charge Q = 410 and quantum number n = 1 (left) and 
n = 2 (right) is shown as a function of the coordinates p = r sin 6 and z = r cos 6. 




FIG. 13: The energy density Tu for fundamental rotating Q-balls with charge Q = 410 and quantum number n = 1 (left) and 
n = 2 (right) is shown as a function of the coordinates p = rsm9 and z = r cos 9. 

rotating Q-balls are exhibited in Table III together with the properties of the corresponding non-rotating Q-ball. The 
properties of radially excited rotating Q-balls have not yet becui studied [4] . 



n 


ids 


J 


M 





0.5927 





293.8 


1 


0.6931 


410 


363.4 


2 


0.7940 


820 


414.7 



TABLE III: Physical properties of fundamental Q-balls with charge Q = 410 and quantum numbers n = — 2. 

We here address the frequency dependence of rotating Q-ball solutions. We focus on fundamental rotating Q-balls 
with quantum number n = 1 and even parity. In Fig. 14 we show the charge Q as a function of the frequency lUs for 
these fundamental rotating Q-balls. We observe the same upper limiting value Wmaxj Eq. (42), for the frequency Wg, 
as for non-rotating Q-balls, which again ensures asymptotically an exponential fall-off of the scalar field cp. For a given 
frequency LUg the charge of a rotating Q-ball is larger than the charge of a non-rotating Q-ball. We thus conjecture, 
that the frequency of rotating Q-balls is also limited by the minimal frequency Wmin, Eq. (43). Unfortunately, for 
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rotating Q-balls the numerical accuracy decreases considerably for large values of the charge Q, refraining us from 
approaching Wmin niore closely. 




0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 2000 4000 6000 8000 10000 

«s Q 

FIG. 14: The charge Q is shown as a function of the frequency u!s (left) and the mass M is shown as a function of the charge 
Q (right) for fundamental rotating Q-balls with quantum number n = 1. For comparison, the charge Q and the mass M also 
shown for fundamental non-rotating Q-balls. The upper branches of the mass M are not discernible (on this scale) from the 
mass of Q free bosons, M — tubQ (right). 

The mass of the rotating Q-balls shows the; same cusp structure as the mass of non- rotating Q-balls. Only the 
minimal charge Qmin of the rotating Q-balls is larger than the minimal charge of the non-rotating Q-balls, Qmin{n = 
1) > Qinin{n = 0). We conclude, that the set of rotating Q-balls exhibits the same general pattern as the set of 
non-rotating Q-balls. 

B. Rotating boson stars 

Wc now turn to rotating boson stars, obtained when gravity is coupled to the rotating Q-ball solutions. We expect 
that the set of rotating boson stars also exhibits the same general pattern as the set of non-rotating boson stars. 
Previously [5, 11] rotating boson stars were obtained (only for a <I>^-potential and) only for a limited range of the 
frequency lUs- Therefore neither a spiral structure nor cusps were observed. 

We here focus on fundamental rotating boson stars with rotational quantum number n = 1 and even parity. We 
exhibit in Fig. 15 the charge Q as a function of the frequency Us for fundamental rotating and non-rotating boson 
stars at gravitational coupling k = 0.2. 

Again, the frequency cOg of the solutions is bounded from above by Wmax, Eq. (42), ensuring an asymptotically 
exponential fall-off of the scalar field. Furthermore, as for the non-rotating boson stars, we observe for the rotating 
boson stars for the smaller values of the frequency lUs a backbending toward larger values of cj^, leading apparently 
also to an inspiralling of the solutions towards a limiting solution. Unfortunately, numerical accuracy does not allow a 
better determination of the spiral and the corresponding limiting values. Comparison of the rotating and non-rotating 
sets of solutions shows that the location of the spiral of the rotating solutions is shifted towards larger values of the 
frequency tOg ■ 

The mass M has an analogous dependence on the frequency lUs as the charge Q, as seen in Fig. 15. When the mass 
M of the rotating boson stars is considered as a function of the charge Q, we observe an analogous cusp structure as 
for the non-rotating boson stars, as illustrated in Fig. 16. For smaller values of k (e.g., k = 0.2) we observe two sets 
of cusps. The first set of cusps is again related to the cusp present in flat space, occurring at the minimal value of 
the charge Qmin, where the upper and lower branch of the rotating Q-ball merge. In curved space the upper branch 
extends again only up to a second cusp, reached at another critical value of the charge, where this upper branch 
merges with a third branch extending down to zero. For larger values of k (e.g., k = 1) this set of two cusps, being a 
remainder of the flat space solutions, is no longer present. 

The second set of cusps is again rcilatcid to the spirals, present only in curved space. Labelling these cusps again 
consecutively A'" = 1, 2, we exhibit the charge Q, the mass M, and the frequency lJs of the boson star solutions at 
the first two cusps in Table IV, for the gravitational coupling constants k = 0.1, 0.2, and 0.3. Numerical problems 
inhibit a better resolution of this cusp structure (with regard to the larger values of N) . 

Addressing finally the K-dependence of the domain of existence of the rotating boson star solutions, we exhibit in 
Fig. 17 the charge Q and the mass M as functions of the frequency Us for several values of the gravitational coupling 
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FIG. 15: The charge Q (left) and the mass M (right) are shown as functions of the frequency uis for fundamental rotating 
(n = 1) and non-rotating boson stars at the gravitational coupling k = 0.2. 
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FIG. 16: The mass M is shown as a function of the charge Q for fundamental rotating (n = 1) boson stars in the full range of 
existence (left), and in the charge range of the first set of cusps (right), for the gravitational coupling constant k = 0.2. 



N 


Qn 


M 


J 




1 


237.2 


179.3 


237.3 


0.4325 


2 


35.2 


66.7 


35.1 


0.9700 


N 


Qn 


M 


J 




1 


111.0 


99.3 


111.0 


0.5100 


2 


23.3 


44.4 


23.3 


0.9700 


TV 


Qn 


M 


J 




1 


910.9 


516.5 


911.1 


0.3400 


2 


76.4 


138.8 


75.7 


0.9450 



TABLE IV: Physical characteristics of the first two cusps of the fundamental rotating (n — 1) boson star solutions for gravita- 
tional coupling constants k = 0.1 (upper), k = 0.2 (second), and k = 0.3 (third). Deviations from the relation J = nQ are due 
to numerical inaccuracy. 



constant k, including the flat space limit. Again, we observe that for values of the frequency Wg close to Wmax, the 
flat space values are not approached monotonically from below with decreasing k. 

We observe from Fig. 17 that also rotating (n = 1) boson stars exist only in an interval [a;o(«;), Wmax]- To obtain 
the lower boimd luq{k) in the limit k ^ oo, we again introduce the scaled scalar fleld 0(r) = ^/K(j}{r), and substitute 
(j){r) in the field equations. Taking the limit k ^ oo, and solving the new set of differential equations, we find for the 
rotating (n = 1) boson stars the lower bound loq{oo) = 0.677. The (scaled) charge Q = kQ of these solutions is shown 
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FIG. 17: The charge Q (left) and the mass M (right) are shown as functions of the frequency uis for the fundamental rotating 
(n = 1) boson stars for several values of the gravitational coupling constant k. Also shown are the limiting flat space values. 

in Fig. 18. Concerning the domain of existence in the limit /t — > 0, we expect a similar pattern as for the non-rotating 
boson stars. For n < 0.06, for instance, wc observe that ojo{k) < Wmin- Numerical problems, however, prevent us from 
obtaining reliable results for very small values of k. 
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FIG. 18: The scaled charge kQ is shown as a function of the frequency uJs for rotating (n ; 
For comparison the scaled charge Q is also shown for non-rotating boson stars. 



1) boson stars in the limit /t — » oo 



VI. CONCLUSIONS 



We have addressed boson stars and their flat space limit, Q-balls. To illustrate the effects of gravity and rotation, 
wc have first recalled the main features of non-rotating Q-balls. Second, we have added gravity to obtain non-rotating 
boson stars. Third, we have added rotation to obtain rotating Q-balls, and finally we have added both gravity and 
rotation to obtain rotating boson stars. Our main emphasis has been to study the general pattern displayed by these 
regular extended objects, and to determine their domain of existence. 

Spherically symmetric Q-balls and boson stars exist only in a limited frequency range. Whereas both mass and 
charge of Q-balls assume a minimal value at a critical frequency, from where they rise monotonically towards both 
smaller and larger frequencies, boson stars show a different type of behaviour. Their mass and charge tend to zero 
when the maximal frequency is approached, while for smaller values of the frequency charge and mass exhibit a 
spiral-like frequency dependence, leading to limiting solutions with finite values of the mass Mum, the charge Qnm 
and the frequency Wiim, depending on the gravitational coupling constant k. 

Denoting the domain of existence of spherically symmetric boson stars by [cuo{K),oJniax], we observe, that when 
K ^ DO, the lower limiting value of the frequency ujo{k) tends to a finite value a;o(oo) < Wmax, obtained by solving 
the set of field equations for the scaled scalar field ^/K(j){r). 

For small values of k, on the other hand, the domain of existence consists of two parts, [wo(k), Wmin] and [o^min, Wmax]- 



17 



In the limit k — > 0, the lower hmiting value of the frequency ll!q{k) tends to a finite value, LUoi^) < <^min- In the interval 
[wo(0), Wmin] there are no flat space solutions, possessing finite charge or mass. The limiting values of the spirals, 
Qiim and Mjim diverge with k^'^/^ when k — > 0. (For a potential with degenerate minima [6] only the single interval 
] needs to be considered, since Wmin = 0.) 

Spherically symmetric Q-balls and boson stars possess radial excitations. Their systematic study has not yet been 
performed. We have, however, considered the first radial excitations of Q-balls and the first radial excitations of 
boson stars at a given gravitational coupling. This indicates, that radially excited Q-balls and boson stars possess an 
analogous frequency dependence and critical behaviour as the corresponding fundamental solutions. 

Rotating Q-balls and boson stars possess a quantized angular momentum, J = nQ. For each rotational quantum 
number n there are even and odd parity solutions. Here we have focussed on solutions with n = 1 and even parity. 
These rotating Q-balls show a frequency dependence and critical behaviour analogous to the non-rotating Q-balls. 
Likewise, the rotating boson stars show a frequency dependence and critical behaviour analogous to the non-rotating 
boson stars. In particular, on the one hand we observe that the mass and charge of rotating boson stars tend to zero 
when the maximal frequency is approached, and on the other hand we see a spiral-like frequency dependence for the 
mass and charge for the smaller values of the frequency. Unfortunately, numerical inaccuracies increase along the 
spirals, preventing a reliable construction of their interior parts. 

Rotating (n = 1) boson stars also exist only in an interval [a;o('*)i (^max]- In the range of k considered, the lower 
limiting value u>q{k) of rotating boson star solutions is larger than the lower limiting value of non-rotating boson 

star solutions. In the limit k — > c», the lower bound ujq{k) is again obtained by solving the set of field equations for 
the scaled scalar field y/K(j){r). 

For small values of k, we also observe an analogous pattern for the rotating boson stars as for the non-rotating 
boson stars. Numerical problems prevent us, however, from obtaining reliable results for the domain of existence for 
very small values of k. 

A systematic study of the frequency dependence of rotating boson stars with higher rotational quantum number 
n and with both parities is still missing and represents a numerical challenge. While the energy density of rotating 
boson stars with n = 1 and even parity has a torus-like structure, the energy density of rotating boson stars with 
n = 1 and odd parity should have a double-torus structure. For radially excited rotating boson stars the structure of 
the energy density may involve concentric tori. Their construction remains an open challenge, as well. 
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APPENDIX A: AUXILIARY FUNCTIONS 



In order to be able to calculate boson star solutions for a given value of the charge Q, we define 2 auxiliary functions, 
p{r,6) and Ws{r,9), The equation of motion for the function p{r,6) consists of the Laplacian acting on p{r,0) and a 
source term, which is proportional to the charge density, Eq. (26), 



P,0 



-l«|V2 I gt^+gtV 



(Al) 



Integration yields 



de + 



dr 



and with 151^2 = l^/'^gr'^ smO/f and Eq. (26) we obtain 

[sinOl^'^r'^p^r) de+ f (sin 6* Z^/^.e) 



dr = 



Q 



Making use of the asymptotic expansion for the function p{r, 9), 



C 



P = Poo + - + 0{ — ]+... , 



1 



(A3) 



(A4) 
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with constants C and poo, and thus asymptotically r'^p,r — > —C, we obtain 

Q 



C / smOde = 2C= , 



This then yields the connection between the function p(r, 6*) at infinity, the frequency ujs and the charge Q, 

= . 



2 Q 

Stt 



(A5) 



(A6) 



Since we cannot impose the value of the frequency Wg, when we impose the value of the charge, we need to solve 
for the frequency ojs- The frequency is a constant. An adequate equation is therefore, that the Laplacian of uJs 
vanishes, 



The appropriate set of boundary conditions for these two auxiliary functions are 



(A7) 

(A8) 
(A9) 



drp\r=0 = , Plr^oo = Poo , dep\e=0 = , d0p\o=^/2 = , 

and Poo is an arbitrary constant. 

When instead of the charge Q a value for the frequency ojg is imposed, use of the auxiliary functions is not 
necessary, but convenient, since the charge is then obtained directly. In this case, the same set of boundary conditions 
is appropriate except for 

(AlO) 



and Woo is the required value. 



APPENDIX B: SYSTEM OF DIFFERENTIAL EQUATIONS 
1. Spherically symmetric solutions 



For the spherically symmetric boson stars we obtain the following set of coupled non-linear ordinary differential 
equations 



d'rf 



1 1 

2 71 



4:Kff UU) - SkloI 0^ + - fWrf - 2ldrf 

r 



+ fdrldrf 
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'2/7 
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-fidrlf 



(Bl) 



(B2) 



d'd = - 



1 1 
2/7 



2u;i 



f 



+ fdrl dr(j) 



For the auxiliary functions the equations are 



= (d^p +^drP+^^ drl drP^ + y 



' d^LOs H drCOs + — drl drijJt 
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(B3) 

(B4) 
(B5) 



2. Axially symmetric solutions 



For the rotating boson stars we obtain the following set of coupled non-linear partial differential equations 
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For the auxiliary functions the equations are 
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APPENDIX C: STRESS-ENERGY-TENSOR 
1. Spherically symmetric solutions 

The nonvanishing components of the stress-energy-tensor for non-rotating boson stars read 




2. Axially symmetric solutions 



nonvanishing components of the stress-energy-tensor for rotating boson stars read 
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